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Effects of non-magnetic disorder on the critical temperature and on diamagnetism of quasi- 
one-dimensional superconductors are reported. The energy of Josephson-coupling between wires 
is considered to be random, which is typical for dirty organic superconductors. We show that 
this randomness destroys phase coherence between wires and that Tc vanishes discontinuously at 
a critical disorder-strength. The parallel and transverse components of the penetration-depth are 
evaluated. They diverge at different critical temperatures Ti^^ and Tc, which correspond to pair- 
breaking and phase-coherence breaking respectively. The interplay between disorder and quantum 
phase fluctuations is shown to result in quantum critical behavior at T = 0, which manifests itself 
as a superconducting-normal metal phase transition of first-order at a critical disorder strength. 
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I. INTRODUCTION 

Although more than a quarter of a century has 
elapsed since the discovery of superconductivity in the 
quasi-one-dimensional (quasi-lD) organic charge-transfer 
(Bechgaard) salts of (TMT^F) aX-type (where TMTSF 
stands for tetramethyltetraselenofulvalinium and X = 
PFq,C104, NO3 being a strong electron acceptor or 
anion) many distinct properties of this material still 
remain a matter of debate. Among those one may men- 
tion the pairing symmetry, the remarkable sensitivity 
of the critical temperature to irradiation'^''*, large up- 
per critical field Hc2 etcetera (see, also Ref.5). The 
irradiation destroys anion order, introducing thus non- 
magnetic damage into system that led Abrikosov to sug- 
gest the possibility of triplet pairing^ in the organic salts. 
Nevertheless, the experimental evidence^ that the spin- 
susceptibility decreases at low magnetic fields in the su- 
perconducting (SC) state of the {TMTSF)2ClOi organic 
conductor, disfavors the triplet pairing mechanism and 
supports spin-singlet pairing. 

Low temperature properties of organic superconduc- 
tors are known to be very sensitive to disorder. Alloy- 
ing with anions, x-ray irradiation, or cooling rate con- 
trolled anion reorientation introduces non-magnetic ran- 
domness into the system, however leaving unchanged, 
to a large extent, the backbone structure and the unit 
cell of the organic superconductors. There is a common 
agreement that disorder, when introduced by means of 
these experimental methods, must be characterized as 
non-magnetic. Yet, it was shown^'^"^" to suppress the 
SC phase. It is worth noting that the methods of gener- 
ating non-magnetic disorder in layered organic supercon- 
ductors are similar (see, e.g. Ref.ll) to those in quasi-lD 
systems, and, therefore a similar mechanism of SC state 
suppression in these two classes of materials may exist. 

Studies of disorder effects on the superconducting 
phase have a long history. The superconducting transi- 



tion temperature Tc for s-wave pairing has been shown to 
be insensitive to the scattering rate on non-magnetic im- 
purities, which is known as the Anderson theorem^^. In 
contrast to non-magnetic impurities, paramagnetic im- 
purities break time-reversal symmetry of the s-pairing, 
and suppress the SC-phase*'^ at some critical concentra- 
tion of the impurities. Correlation between the para- 
magnetic impurities via the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction yields a spin-glass phase be- 
low T — Tg, which was shown^^ to shift the supercon- 
ducting transition point towards higher temperatures. 
The Anderson theorem is not applicable to unconven- 
tional superconductors with d-wave pairing symmetry 
A small concentration of non-magnetic impurities may 
destroy d-wave pairing, producing a finite lifetime for 
quasi-particles near the nodes in the gap^^^^^. 

Inter-chain (inter-layer) pairings as well as intra-chain 
(intra-layer) pairings, occurring at different local critical 
temperatures in quasi-lD (quasi-2D) systems, yield also 
an inhomogeneous nodal order parameter^^'^", which af- 
fect considerably the upper critical magnetic field. 

Suppression of superconductivity in the presence of 
non-magnetic impurities can in general be realized by 
destroying either the modulus or the phase coherence of 
the order parameter. Interplay between superconductiv- 
ity and Anderson localization in a strongly disordered 
superconductor was shown^^"'^" to result in spatial inho- 
mogeneity of the order parameter. Diffusive scattering 
of particles in the random field of non-magnetic impuri- 
ties enhances Coulomb repulsion^^'^^, and consequently, 
reduces the amplitude of the order parameter. Meso- 
scopic fiuctuations in a superconducting thin film were 
also shown'^^ to yield a spatial inhomogeneity of the or- 
der parameter. 

Effects of order parameter phase fiuctuations on the su- 
perconducting transition temperature Tc have been stud- 
ied in Refs. 34-38. It is well known that there is no 
SC phase transition in one-dimensional (ID) and two- 
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dimensional (2D) systems"^^, since strong fluctuations 
of the order parameter phase destroy off-diagonal long- 
range order (ODLRO) in a single SC wire and in an iso- 
lated SC film. Even a small interchain-coupling in clean 
quasi- ID superconductors restores however ODLRO to- 
gether with a finite transition temperature. The suppres- 
sion of Tc by strong phase fluctuations in clean quasi- 
ID superconductors was analyzed in Ref.34,35. Classi- 
fying superconductors with small stiffness as bad met- 
als, Emery and Kivelson have evaluated^'' a critical tem- 
perature of phase ordering by formally dividing a clean 
bulk superconductor into small regions with well defined 
phase, and have shown considerable suppression of an 
SC-phase by phase fluctuations. Effects of disorder on 
phase fluctuations were however not considered in all of 
these papers. 

Distinct structural peculiarities of quasi- ID organic su- 
perconductors demand a special theory, which should 
take into account non-magnetic randomness as well as 
phase fluctuations. Indeed, the high purity of the organic 
superconductor backbone, even in an overall dirty limit, 
excludes spatial inhomogeneity of the order parameter 
modulus along the SC-wires. This renders inapplicable 
the above-mentioned theories for pair breaking. 

In contrast to these previous activities we study 
in this article a disorder-driven superconductor-normal 
metal phase transitions due to the competition of non- 
magnetic randomness and phase-fluctuations in quasi- ID 
superconductors. We consider weakly linked quasi-one- 
dimensional superconductors with random Josephson- 
couplings between pure one-dimensional (ID) SC wires. 
Singlet pairing is assumed within the wires. Therefore, 
we assume that non- magnetic randomness does not affect 
the order parameter amplitude. 

We demonstrate in this article that 

(i) non-magnetic randomness in the Josephson- 
coupling destroys correlation of the phases between dif- 
ferent chains in quasi- ID superconductors even in the 
classical phase fluctuation regime, 

(ii) non-magnetic randomness yields quantum critical 
behavior in addition. A superconducting normal-metal 
phase transition occurs at T = increasing the strength 
of disorder, and 

(iii) a suppression of the SC phase occurs discontinu- 
ously as well both the classical and the quantum phase 
fluctuation regimes as a first-order phase transition when 
the disorder-strength reaches a critical value. We derive 
parallel and perpendicular components of the penetra- 
tion depth, A|| and A_l, and show that they diverge at 

different critical temperatures tJ^' and Tc, which corre- 
spond to pair-breaking in the wires and to phase coher- 
ence breaking between the SC wires, respectively. 

The paper is organized as follows: In Sec. II we study 
the interplay of randomness in the Joscphson energy 
with phase fluctuations inside the classical regime. Al- 
though an arbitrary small interchain coupling in a clean 
quasi-lD superconductor stabilizes the ODLRO giving 
a finite transition temperature Tg, we show in this Sec- 



tion that non-magnetic disorder in the Joscphson energy 
suppresses discontinuously when the disorder-strength 
reaches a critical value. In Sec. Ill, we investigate effects 
of randomness on the transverse rigidity and on in 
the quantum fluctuations regime. We show in this Sec- 
tion that a suppression of the SC phase is managed by 
two parameters characterizing the disorder-strength and 
a dynamical charging parameter a in the system. The 
quantum criticality at T = and the phase transition at 
flnite temperature are described separately. The trans- 
verse rigidity in the field of the phases is shown to vanish 
discontinuously. The jump at the breakdown point de- 
creases monotonically with increasing a in the interval 
of < a < 1, and vanishes for a = 1. The breakdown 
point is pushed towards higher values of the disorder- 
strength as a increases. We show that, under a change 
of the disorder-strength, the critical temperature evolves 
similarly as the transverse rigidity at T = 0. The anal- 
ysis of diamagnetism in this Section reveals completely 
different behavior for parallel and perpendicular compo- 
nents of the penetration depth. The penetration depth 
for a magnetic field parallel to the SC-wires is shown 
to be defined by the phase-phase correlator between two 
neighboring wires, which behaves non-linearly in tem- 
perature and reveals a discontinuous dependence on dis- 
order; by contrast the perpendicular penetration depth 
does not depend on disorder and shows a conventional 
temperature dependence. In Sec. V we summarize re- 
sults obtained in the paper. The explicit evaluation of 
the phase-phase correlator is given in the Appendix. 



II. CLASSICAL PHASE FLUCTUATION 
REGIME 

A quasi-lD superconductor is modeled as a system of 
one-dimensional wires, which are placed regularly and 
parallel to each other, forming for example a square lat- 
tice in the cross-section. Weak tunneling between the 
chains results in an open Fermi surface for the normal 
metallic state, and yields also the Josephson-coupling be- 
tween nearest-neighboring chains in the superconducting 
state and strong anisotropy in kinetic properties. 

The free energy functional of a quasi- ID superconduc- 
tor, which is weakly linked with Josephson-coupling en- 
ergy -Ej,j-i-g between nearest-neighbor chains, can be writ- 
ten in the presence of the magnetic field B as 



2ee|| 



Smii^M \ dz he 
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where <pj(z) denotes the phase of the order parameter 
Aj(2;) = |Aj| exp(i(pj(z)) at a point with dimensionless 



3 



coordinates r = {j = {jx,jy},z}, A = {Aj_,A^} is the 
vector-potential, and A^i^-* the linear density of SC elec- 
trons with nP{T) = 7vi^^(0) = N^^^ ~ at T > T^, 



at T < 



and N^^\t) = M^\0)t{T) with t{T) = 

Tc^\ Dimensionless coordinates r = {j, z} are introduced 

on the scale of the longitudinal = ^ and the 

transverse ^ a_L components of the coherence length, 
so r = {ry = {jxa±,jyajL},z} -i> {j = {jx,ij,}, 
The last term in Eq.(l) describes the Josephson-coupling 
with the coupling energy i?j,j+g between the wires, which 
is minimal for a coherent tunneling {(pj{z) ~ (pj-^-g{z)) of 
Cooper pairs between neighboring wires. 

Fluctuations of the order parameter modulus can be 

neglected for pure superconductors'^^ far from Tc^\ sat- 
isfying the condition (ri^' - T)/T^^^ > n-^/s^ where 

Ti^^ is the mean-field critical temperature calculated for 
an isolated wire and n is the number of bands in each 
chain^*^. Therefore, the contributions to the free energy 
functional (1), coming from the modulus of the order pa- 
rameter |Aj|, are omitted. 

We assume the Josephson energy -Ejj+g to be a random 
parameter with Gaussian distribution, centered at the 
mean value E^, given by 
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The variance W"^ is taken as a measure of disorder 

strength in this coupling of nearest neighbor chains. Em- 
ploying the replica trick one can calculate the average 
value of the free energy = — T(lnZ) over disorder. 
As usual we use (InZ) = lim„^o ^(^") ^.nd, in addi- 
tion, express the n-th power of the partition function 
Z = / rij ^"^jS"^''^^ by means of replicated fields tp", 
a = l...n, as (for B = 0) 



in 



+ 5^Eg[l-cos(v?-^]Vg)]l+^E 

_ o 



8to||^^ V dz 



The quadratic cosine term is linearized with the help of 
a Hubbard-Stratonovich transformation by introducing 

an auxiliary field g. As a result, the sum over the 
replica variable a is factorized and the replica limit can 



be performed, yielding for the averaged free energy 



with F 



j,g 

xlny Uv^^e-^/^ 



2 ""J.g X 



(4) 



+ 



E(^s - N('k\\WCi,s)[l - cos(^j - <^j+g)] j. 



(5) 



The average value of a given functional C({<^j}), e.g. 
cos (pj or cos((^j — <^j+g), can be obtained according to the 
relation {{C({^i}))) = -r^(lnZ)|^.=o by adding the 
source term J2j I dzij^C {{(fi^}) to the free energy func- 
tional, which yields for the correlator 



{{CiWi}))) 



-n- N^s^kll 



2tt 



(6) 



where the double bracket ((...)) means averaging over 
thermodynamic fluctuations and over randomness. In or- 
der to estimate an asymptotic behavior of the correlator, 
e.g. ((cos(i^j — v'j+g))) we write the integrand of Eq.(6) 
as exp{— Afi^^^^|/(Cj^g)}, and apply the stationary-phase 
approximation to determine an extremal vahie of the aux- 
iliary field Cj^g minimizing the function /(Cj,g)- The min- 
imal value of Cj,g is obtained to be 







J d2;|(cos(<pj(2;)-(pj+j 



(cos(yj(2:) - (Pj+sjz)) cos(yj(0) - yj+g(0))) \ , . 

(cos((^j(0)-¥'j+g(0))) /■ ^' 



The constant on the exponent can be estimated to 

be equal to ivi"'^''^|| ~ ^ 10"^ for the organic supercon- 

2 ductors with being the Fermi energy, which ensures a 
sharply peaked saddle point of the integrand. The ther- 
modynamic averages in the expression of Cj.g are taken 
with the free energy functional, given by Eq.(5), at the 
saddle point Cj,g = Cj,g- So, a contribution of the non- 
magnetic randomness to the effective free-energy func- 
tional is proportional to the variance of the phase cor- 

(3)relator, which gives an idea on the form of the disorder- 
dependent term in the effective functional. Note also that 
the saddle point for the averaged free energy is given 
as Cj,g = / dz(^^ ^ cos(^j(z) - ^j+g(z))]), where J" 
is the value of the free energy at the saddle point. 

The critical temperature for the quasi- ID supercon- 
ductors can now be found from Eq. (6), written for cos(/?j 
by using the self-consistent mean-field method^**, which 
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consists in replacing the phase correlations of the cosine 
term by 



i?^[l-((cosM))cos(^(z))], 



(8) 



where E± — Eg-^g- ^ clean system ((cos((/3)))e// 
was chosen^'' to be equal to (cos((^)). For the disordered 
superconductor we choose {{cos (p))eff = {{cos(p)) — 
^a)^^^ ((cosy)V((cosy))'' ^ The functional integral over 

the phases in Eq.(6) can not yet be evaluated, even after 
this simplification. Taking advantage of the smallness of 
{Ej_ — Ns^^ ^\\WQ{{cos{ip))) ef f near T^, we expand both 
the numerator and the denominator of the integrand of 
Eq.(6), written for ((cos((^)))e//, in this parameter. The 
thermodynamic averages become pure one-dimensional 
after this expansion, which can be taken easily, yielding 
a power series of C for the integrand. Therefore, the inte- 
gration over C. is immediately performed. Since all higher 
order in ((cos(iy9)))e// terms of the expansion vanish at 
T = Tc, we get the equation for Tc 



X {cos{ip{0)) cos{ip{z)))dz, 



(9) 



where r] is the coordination number. The phase-phase 
correlator in Eq.(9) is calculated in the clean limit of 
the ID free energy functional, obtained from Eq.(l) by 
setting i?j j+g — 0, which returns (see for example Ref.^^) 

(cos((^(0)) cos{ip{z))) = exp{-|^|/re}, (10) 

where = fi^NP {T)/2m\\^\\T. Introducing a dimen- 
sionless T^-shift by 



t = ^/r^epE^ ( Y 



1 



(1) 



(11) 



with e F being the Fermi energy, and a dimensionless dis- 
order parameter 




Eqs.(9), (10) yield 
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1 = t^{l - qt). 



(1) 



Ei_ ' 



(12) 



(13) 



The full solution of Eq.(13) has three roots, among 
which the physical one is confined to the finite g-range 
as shown by the bold line in Fig.l. A (physical) solution 
exists thus only within the finite range between clean 
limit (CL) and dirty limit (DL). One may expand and 
control this physical solution (of the cubic equation) in 
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FIG. 1: The physical solution t{q), giving the Tc{W)- 
dependence, within the full range from clean limit (CL: g = 0) 
to the dirty limit (DL: qc = 2/3-^/3)) is highlighted as the bold 
(blue) curve. Formal solutions of the cubic Eq.(13) are shown 
for completeness. Tc{q) vanishes abruptly at q — qc. 



the weak disorder regime (small q), where the Tc-shift 
obeys 



1 

TV. 
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showing that Tc decreases with increasing randomness 
like . For a pure system Eq.(14) gives the dependence 

1/2 

Tc ^ , in agreement with Efetov and Larkin in Ref. 
34. This expression shows that even a small interchain- 
coupling sustains ODLRO in the system and, conse- 
quently, the critical temperature increases with y/E±. 
On the other hand, the competing destructive effect of 
disorder reduces Tc due to " melting" of the order param- 
eter phase coherence between neighboring chains. 

According to the (physical) solution of Eq.(13), the 
critical temperature decreases monotonically with in- 
creasing q but finite Tc are confined to the interval 
< g < 9c = The SC-phase becomes fully sup- 

pressed for disorder-strengths W'^ exceeding a critical 
disorder- value Wc given by 



^E^Tc 
3 



(1) 



lEj_ 



(15) 



beyond which the system is in a normal metallic phase 
(for W'^ > Wc)- The critical temperature drops to zero 
at — Wc with a jump of size 



ATc = T* = 



rjepEj 



1 



T 



(1) 



(16) 



Thus the SC-normal metal phase transition appears as a 
Ist-order transition. 

In order to describe the behavior of t near the disor- 
der limit DL (see Fig.l) we expand around {t*,q*} = 
{VS, in terms of small (nonnegative) St = t* — t 

and Sq = q* — q, which gives 6t = 3^^^\/Jq. 
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Near the dirty limit, the Tc-variation has an infinite 
slope (see also Fig.l). This can be reexpressed in terms 
of the physical parameters ST^ = Tc — T* and the 
variance of the Josephson-coupling, by reinserting 
Eqs.(12,14), as 



in the vicinity of the breakdown point 



{W^ - W^f/^ (17) 



+ 



-.(1) 



(1) 



Thus the critical temperature decreases with disorder al- 
most linearly but, approaching the dirty limit DL, it 
finally turns into a (nonanalytic) square-root behavior 
close to the breakdown point. In the absence of the 
non-magnetic disorder, even arbitrarily small Josephson- 
coupling between the chains stabilizes the ODLRO and 
gives a non-zero critical temperature. However, the SC 
phase with finite Josephson-coupling can be destroyed by 
increasing the strength of non-magnetic disorder. 



III. QUANTUM PHASE FLUCTUATIONS 

The self-consistent mean-field method, applied above 
for the classical phase-fluctuation regime, expressed the 
Tc-equation in terms of the pure ID phase-correlator 
(10), neglecting in this respect the Josephson-coupling 
between neighboring chains. We shall now improve the 
calculation of the phase-correlator by taking into account 
the transverse rigidity of the system, which provides a 
more realistic determination of the transition tempera- 
ture in the quantum fluctuation regime. Our calcula- 
tions are carried out in the Hamilton formalism for con- 
venience, yet the problem can be formulated in the path 
integral language^^ as well. 

Let us start from the Lagrangian, again taking B = 
for simplicity, 



ifeii(o) 



Kt"^}, (18) 



where (p denotes the time derivative of the phase. The 

dynamical term in the Lagrangian can be interpreted as 
the electrostatic energy of charged wires^^'^® 

= ^ E / / - ^'mmiz'), (19) 



generated according to the first Josephson equation (p = 
{2e/h)V, and Cij(2; — z') are the specific coeSicients of 
electrostatic induction. Rewriting the electrostatic en- 
ergy Ef,i in terms of the time-derivative of phases, the 
Fourier transform /^(qj,,^^) of the coeSicients Ki^j{z — 



z') = -^Cij{z — z'), has the physical meaning of a com- 
pressibility. In Eq.(18) we neglect dispersion in the com- 
pressibility and assume K{q±,qz) = K = const. This 
approximation is equivalent to a replacement of the elec- 
trostatic energy (19) by / dzJ2j ^^i^)- The param- 
eter K can be calculated^^ in the presence of Coulomb 

2 

screening for a small Born parameter < 1, which 

results in 



K 



nhvo 



TTfivo a 



n -1 



(20) 



There, Kq = is the unscreened compressibility, Vq 
denotes the longitudinal velocity of an electron on the 
Fermi surface averaged over n siibbands, while aj_ and d 
stand for the interchain-distance and the diameter of a 
superconducting wire, respectively. 

In order to get F!:/-^{ip} we average the free energy 
over disorder and apply a mean-field approximation, cor- 
responding to the replacement J dz J dz'[l — cos{(fj{z) — 
'P}+s(zM^-cos{ipj{z')-(fij+g{z'))] by (([l-cos(<^j(2;)- 
"Pi+siz))])) I <iz[^ - cos(¥'j(2) - ^j+s{z))]. 

The Hamiltonian, expressed through the phases and 
canonical conjugate momenta Hj , becomes 



n = Y:ldz[ 



dz 



E'^cJl - cos{(pj{z) - ipj+g{z))] 



+ 



(21) 



where Uj = j-j^ = ifi-i^^y (0)(/jj, while 6ci is given by 



E±T 



{{[1 - COs{ipj(z) - ipj+g{z))])) 



(22) 

and reprc;sc;nts the diniensionless anisotropy-parameter 
or the transverse rigidity of the system; So in Eq.(22) is 
transverse rigidity of the pure system 



^0 



Ei_ 



1/2 



{epEA,) 



1/2 



■.(1) 



(23) 



The phase dynamics in the classical limit can be ob- 
tained from the Euler-Lagrange equation, which is de- 
scribed by a set of coupled sine-Gordon type non-linear 
equations 



^j(^) 



dz^ 



^52sin((pj(z)-^j+g(z))j, (24) 



where a; is a characteristic scale of frequency, and 



iVf^%) 



(25) 
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We express u) as oj = 2TraT^/'^TP/h, where 



a 



16 



1/2 



(26) 



The parameter a is the essential parameter of the the- 
ory, which can be written, using Eq.(20), as 



a 



1 + 



-,1/2 



(27) 



The system of equations (24) is hnearized for small ipj 
and its Fourier-transformation becomes diagonal 



<fi{fi±,qz) = -w^(q_L,5^)v'(q-L,fe) 



(28) 



The eigenfrequency of oscillations u{qj_,qz) is given in 
the harmonic approximation as 

io{ci, q,) =u)[q''^+Sli2{2- cos q^- cos qy)]^/\ (29) 

These equations describe the frequency of low-lying 

plasmon-modc of the system. 

The quantum description is realized by expressing ipq 
and Ilq as a linear superposition of Bose operators 6q and 
bl [6q,6^] = l,by 



w(q) 



1/2 



1/2 



If we expand the cosine term in Eq.(21) up to the 
quadratic term and express the phases and the conju- 
gate momentum Ilq in terms of creation and annihilation 
operators, we get the Hamiltonian in the harmonic ap- 
proximation as 



Ho 



q 



(31) 



where the energy spectrum is defined by Eq.(29). 

In order to take into account the quantum effects in 

the Hamiltonian, we have to express the cosine term in 
Eq.(21) in a normal ordering before expanding over 
and 6a 



1 - 
where 



l-C0s(<Pj(2)-(^j+g(z)) = 

2 V 



h.c 



(TTOUJ 



1/2 



(1 



(32) 



The prefactor cxp{— 5a(g, 0)} originates in the com- 
mutation relation between 6]^ and bq, and by taking into 



account the Baker-Campbell-HausdorfF relation exp(J?i + 
H2) = exp(ffi) exp(ff2) exp{— ^[iJi, 1/2]}- Furthermore 



Sa{g,0) 



-E 

2 ^ 
q 



^q(j,g)l' = 



naoj ^ 1 — cos(q_Lg) 



N 



E 



w(q_L,52) 



(34) 

where A'' is the number of unit cells per volume. 

It is clear from Eq.(32) that the physical meaning of 
exp{— ^^(g, 0)} is an average of cos(</?j(2) — <^j-|-g(2;)) over 
the ground state at T = 0. In the framework of the self- 
consistent plionon approximation (SCPA), we expand the 
Josephson term (32) in powers of the creation and anni- 
hilation operators, 6^ and 6q repsectively. Expressing 
the leading (harmonic) part of the Hamiltonian (21) in 
terms of the particle number operator Nq = bl^bq, we 
obtain again a harmonic Hamiltonian as in Eq. (31), 

?io lia^O), but with the different oscillation fre- 
quency 

a;(q, q.) = w[ql + ,5,^-^" '^^'°^2(2 - cosg, - cosg^)]^^ 

(35) 

Thus the application of the SCPA results in a re- 
normalization of the parameter of anisotropy dd in the 
frequency of oscillation a;(q, 0) by means of the phase- 
phase correlator as 



(36) 



We note that in order to preserve a maximal coher- 
ence of phases at ipj = fj+g in the Josephson term 
l—cos{(fj{z) — (f>j+g{z)) the latter is replaced in the frame- 
work of the SCPA by e-^°(s,o) _ cos((^j(^) - 0j+g{z)), 
which corresponds to shifting of the energy origin. In- 
deed, the Josephson term in the initial expression of 
the Hamilton function was introduced in a such way 
that it becomes zero for a maximal coherence of phases 
(pj = (fj+g- Zero-point fluctuations at T = in the quan- 
tum case destroy the phase coherence and increase the 
Josephson energy. By shifting 1 — )■ e"'^"^^'^) we again 
reach a minimal Josephson energy at T = in the quan- 
tum case too. Similar shifting was done also in Eq.(22). 

In the expansion of the exponential operator of Eq. (32) 
we select all diagonal terms, which can be expressed in 
terms of the bosonic particle number operator Nq = 
6^6q38. This yields 



P-Sc(g,0) 



n - 



\A^\XK + lAq^l^^ - l^ql^^^ ± . . . ^(37) 



2! 2! 



3! 3! 



Absorbing now the product in Eq.(37) into the exponen- 
tial form, and neglecting all higher orders in |Aqp beyond 
the leading term |AqP (for justification see below [41]), 
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we get 



where the quantum parameter of randomness Qgu reads 



{C0s((^j - <Pj+g)}diag = exp{-5„(g, 0) - ^ l^ql^-^q} = 



e-Sa(g,T)=e^p. 



2TTau 1 — cos(q_Lg) 



N 



■E 



w(q) 




Thus, after this step of calculations we still restrict our- 
self with harmonic approximation, describing the system 
by the Hamiltonian 'Ha\T) like (31) where the trans- 
verse rigidity d'^i in the frequency a;(q, T) is renormalized 
as 

Sli ^ Sl^{T) = 5l exp{-S„(g, T)}, (38) 



where 



-S„(g,T) 



((C0s((pj(z) - (^j+g(z)))) 



Tr{e-^^^ 'cos(yj(^)-yj+g(^))} 



Tr{e 



(0) 



} 



(39) 



The trace over the diagonal part of the phase-phase 
correlator within the harmonic approximation replaces 
the bosonic filling number operator TVq by the Planck dis- 
tribution function for phonons with energies of huj{ci, T) 
as TVq {exp (^ ^(q-'^) j _ yielding the following 

expression for the 8^ (g, T) 



nau) 



E l - eos(q_Lg) f fiu{ci, T) 



(40) 

The correlator e •^"(g.r) y _ 0-limit, as given 

by Eqs.(40) and (34) respectively, are evaluated explicitly 
in the Appendix. 



A. Quantum Criticality at T=0 

The zero-temperature behavior of the system is ana- 
lyzed by means of the phase-phase correlator e~^°(S'"\ 
the explicit expression for which is given by Eq.(59) 
in the Appendix. Expressing the phase-phase corre- 
lator e~'^°(S'°) in terms of Sa{g,0) gives e~'^°(8'°^ = 
(^gM(O))" = S"^, which implies that even a small 
interchain-coupling stabilizes ODLRO in the system, 
hence also a finite T phase transition should exist. In 
order to get an explicit expression for the dependence of 
Squ on So and on disorder, we have to solve the equa- 
tion 5^^ = <52^e-'5°(8,o) together with Eq.(22) for Sd- 
Thus the equation for the reduced transverse rigidity 



^qu ~ ^qu/^'qu, whcre Sq'„' = 5^^" \s the renormalized 
transverse rigidity for the clean system at T = 0, as- 
sumes the form 



r(0) 



Qguj 



(41) 



Qqu 



(42) 



C in Eq.(42) is a constant C ~ 1. Although Eqs.(13) and 
(41) are written for two different characteristic parame- 
ters of the system, it is easy to see that the equation 
for S*^{0), if we neglect the quantum effects at a = 0, 
coincides with Eq.(13) written for y = 1/t. 

Eq.(41) can be solved approximately for moderately 
weak disorder, yielding the following expression for Sgu (0) 



Squ{0) 



S72 



1-C- 



w 



So- 



j/2 



(43) 

Hence, the evaluation of the phase-phase correlator 
{{cos{ipj{z) — ^j+g{z))))\T=o = e""^"'^'*^^ in the presence 
of disorder yields the result 



{{C0s{ipi{z) - ipi+g{z))))\T=0 



( epSA 



(44) 



In the absence of the disorder, i.e. for W = 0, we retrieve 
the phase-phase correlator 



_ ( epEA 



as obtained by Efetov and Larkin^"* for pure quasi-lD su- 
perconductors. The critical temperature for a quasi-lD 
superconductor, according to Efetov and Larkin^^, is de- 



fined by Teo = 4«(0)Tc^'' for pure superconductors and 
by Tco = 5qu{Q)Ti^^ in our case for dirty superconduc- 
tors. This relation with Eq.(44) shows that Tco decreases 
nonlinearly with disorder. 

The numerical solution of Eq.(41) is depicted in Fig. 2. 
The reduced T = transverse rigidity <5q„(gq«)|T=o is 
shown to decrease with increasing disorder for (fixed) 
a < 1, and suddenly drops to zero at the critical disorder 
strength = Hence the quantum critical behav- 
ior corresponds to a first order phase transition. Fig. 2 
shows how the breakdown point shifts with increasing a 
to higher values of the randomness, and the jump van- 
ishes as a — > 1. Eq.(41) becomes linear for a = 1 and 
gives, by inferring the (M^)-relation from Eq.(42), 



Squ{W)\T=0,a=l = 5^ 



1 



CW^ep 



(45) 



Here, the transverse rigidity 5qu{W)\T=Q,a=i decreases 




0.00 0.05 0. ! 0. 1 5 0.20 0.25 0.30 0.35 



qqu 



0.10 0.15 0.20 0.25 

Residual Resistivity (H em) 



FIG. 2: The dependence of the reduced T — transverse 
rigidity Sg^{qqu) on the disorder-strength parameter Qqu is 
shown for < a < 1 in steps of Aa — 0.1. At qqu = Qqn, 
5qu{<lqu) drops to zero for a < 1 and vanishes continuously 
only at a = 1. Inserts show the a-variation of the jump 
(upper right corner) and of its position gq„(a) (lower left). 



linearly with increasing and vanishes at 




The quantum critical behavior in the model is however 
controlled by two parameters, the strength of randomness 
QquiW) and the parameter of quantum dynamics a. For 
a < 1, the superconductor- normal metal phase transition 
at r = is always discontinuous, and only turns into 
second-order at a = 1. 



B. Phase transition at finite temperatures 

Let us now study the finite T behavior of the trans- 
verse rigidity. The phase transition in a quasi- ID su- 
perconductor occurs at some temperature T — Tc when 
the transverse rigidity in the ensemble of phases {Vj(^)} 
vanishes, which results in "melting" of the phase coher- 
ence. The energy spectrum w(q_L, qz) of the collective ex- 
citations is reorganized and the transverse q_L-dependent 
part of CL'(q_L,qz) vanishes at T = Tc, i.e. symmetry 
breaking occurs in the bosonic excitation at T = Tc. 
Inserting the solution of Eq.(34) for T < aTi^-* into 
5l^{T) = ^2 g-s<,(g,T) j^^j ^gjj^g Eq.(38), we obtain 

^S.(r.^.?„(0)(^)°cxp{-cXMi},,„, 

where a new temperature scale is introduced by means 
of Tco = <5,„(0)Ti^\ In terms of y = {^Y'^ MO 

and = (^) ^ Eq.(47) assumes the form 

y = exp{—6/y}, which has a non-zero solution only for 



FIG. 3: Comparison of the Tc-decrease, as obtained from 
Eq.(48) with fit parameter a = 0.47 of the present theory (af- 
ter proper rescaling), with experimental data imported from 
Ref.lO by Joo et al. Different symbols belong to different 
samples and different slow (squares) or fast (triangles) cool- 
ing procedures. In a first approximation the theoretical solid 
curve, as shown from the clean Umit up to the dirty limit, fits 
and confirms the slightly but increasingly non-linear behavior 
with a final discontinuous drop of Tc at the critical disor- 
der. In the lower left inset Tc-curve employing the Abrikosov- 
Gor'kov-Larkin (AGL) digamma-function (dashed curve) for 
unconventional superconductors is shown for comparison with 
the present theory: the AGL-curve is chosen such that its 
value and slope at zero disorder agrees with the present one. 
The insert in the upper right corner shows the remarkable 
deviation of these two theoretical curves in the large disorder 
regime and close to the breakdown. 

9 < e^^ . The finite solution of this equation vanishes dis- 
continuously ai 9 = 6c ~ e^^, giving the following value 
for Tc 

rc = <59n(0)TWa-^(2/eC)^. (48) 

The magnitude of the jump in y{9c) is e~^, and hence 
the phase transition is of first-order. The similar behavior 
has been found also in the planar rotor model^^ in the 
absence of disorder. The dependence of on the dis- 
order is determined by the zero-temperature transverse 
rigidity (5qu(0), the behavior of which is depicted in Fig. 
2. Therefore, for arbitrary a < 1 the critical temper- 
ature decreases monotonically with increasing random- 
ness and drops to zero at the critical disorder strength 
<Iqu — Qqu- The variation of the critical temperature ver- 
sus the residual resistivity of the organic superconduc- 
tor {TMTSF)2{Cl04)(i-^){Re04)x has been explored 
by Joo et al. in Ref.9,10. The experimental data are 
read off from Ref.lO and copied into Fig. 3 in order to 
provide a close comparison with a theoretical fit-curve as 
obtained from Eq.(48) of the present theory. 

In this first approach and within a moderate accuracy, 
the data as published in the paper of Joo et al., eventually 
appear to find an explanation by our theory. In order 
to establish the link between experiment and theory, the 
following argument is exploited: the substitution of CIO4 
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anions by ReO^ in the relaxed (R) samples increases the 
residual resistivity, which is proportional to the inverse 
lifetime 1 /r of electron as well as to the disorder strength 
W"^ (as used in the present theory). 

For the relaxed samples of Ref.lO, referred to as the 
R-sample(s), the doping concentration x varies in the 
interval of < x < 0.1 under slow cooling from the 
clean sample to the nominal concentration x = 0.1. Slow 
cooling of i?-samples assures a uniform orientation of 
the anions along the stacking axis, whereas fast cooling 
in the (quenched) Q-samples introduces strong orienta- 
tional disorder and increases the residual resistivity. Tc 
decreases quasi-linearly with increasing disorder (or the 
residual resistivity) in the large interval of the random- 
ness. Around the breakdown point the dependence of Tc 
on randomness is non-linear. The critical doping concen- 
tration, corresponding to the breakdown of the supercon- 
ducting state, grows with the quantum charging effect in 
the system. Substitution of CIO4 anions by Re04 seems 
to increase the quantum charging parameter a, shifting 
thus the breakdown point to a higher value of the residual 
resistivity. All these features and experimental evidences 
agree well with the theory. 

In order to compare with other well-known cases of 
Tc-suppression by disorder we consider the pair break- 
ing theory for a superconductor with unconventional 
gap symmetry^^'^*. This physically different case of Tc- 
reduction by non-magnetic impurities in unconventional 
superconductors was found to be described by the fa- 
mous digamma-formula of Abrikosov-Gor'kov's conven- 
tional pair breaking theory^'^ in the presence of paramag- 
netic impurities. The Tg-reduction may thus be expressed 
in the form 



In 



= * 




(49) 



where ^ means the digamma function, p = h/2TTc^'' 
is the depairing parameter, and r the elastic scattering 
time. Both theoretical curves, as shown in the insets of 
Fig. 3, can be chosen to coincide for weak and moderately 
strong disorder (where the linear decay is rather unspe- 
cific) . Approaching the SC- breakdown at larger disorder 
they differ however substantially. The AGL-solution for 
unconventional pairing approaches Tc = continuously 
and obeys a square root dependence Tc ^ (q — 9c) 
where q stands for the disorder-strength in the AGL-case. 
This square root law follows from the leading (0(z^)) cor- 
rection of the digamma- function \l/(i -|- |) given by the 
Laurent series of its exponential 



exp 



- + ^-+0{z^) (50) 
z 24q ^ ^ ^ ^ 



near the logarithmic branch point of ^ at infinity, hence 
z = 0. In the physical context the variable z corresponds 
to the critical temperature Tc of Eq.(49). By comparing 
the exponential of Eq. (49) one can see that z = Tc = 

is reached for q = pT^^^ /{2i:) qc = exp(*(i)). 



By contrast, the present theory does not allow for a 
continuous breakdown of superconductivity. According 
to Fig. 3 the suppression of Tc is stronger and an abrupt 
breakdown occurs at Tc\min > 0. The numerical data 
show a square root behavior however near the minimal 
finite Tc- 

The breakdown point in the curve, corresponding to 
our theory, seems to allow for the existence of an inter- 
mediate phase, perhaps a glassy phase below a tricritical 
point. 



IV. MEISSNER EFFECT 



The current density is calculated according to 
ij(z,j) = -T^{lnZ{A)), where Z = /Pc^e"^-/^. 
The complete expression for J in the linear response ap- 
proximation can be obtained after averaging of In Z over 
disorder (2) by using Eq.(l) for Fgt- One obtains 



^11 



2mii dz 



m\\c 



(51) 



for the longitudinal component of the current and 



2ea_L^I|A^i 



(1) 



E t^7^S^«((cos(¥'j - ^j+g)))(gA^) (52) 



for the transverse component of the current. 

For simplicity we present here only the diamagnetic 
contribution to the i-th {i = ||, -L) component of the cur- 
rent 



(53) 



where the longitudinal (Ay) and the transverse compo- 
nent (A_l) of the penetration depth are obtained as 



(54) 



and 



2 &TTe^Ni''\T)E^ 



((cos((^j-<Pj+g))). (55) 



While A||(r) diverges at T = t]^^ due to pair break- 
ing in the SC wires, A_l(T) diverges at the global SC 
transition temperature T = Tc, where the phase coher- 
ence between neighboring wires is destroyed. The tem- 
perature and the randomness dependencies of A|| (T) and 
A_l(T) also strongly differ each other. The transverse 
component of the penetration depth is determined by the 
phase-phase correlator, revealing non-linear temperature 
dependence and discontinuous behavior at the critical 
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disorder strength. Nevertheless the longitudinal compo- 
nent of the penetration depth is given by the conventional 
London expression and does not depend on the disorder 
strength. Randomness in the Josephson coupling shifts 
Tc to lower temperatures and, therefore, the magnetic 
field parallel to the SC wires penetrates easier into the 
organic superconductor. On the other hand, the type 
of disorder considered in this article does not break the 
Cooper pairs, keeping thus the penetration of a perpen- 
dicular magnetic field into the SC wires unchanged. 



V. CONCLUSIONS 

In this paper we report disorder-effects on Tc and on 
the diamagnetism of quasi- ID superconductors with ran- 
dom Josephson-couplings. Interplay of non-magnetic dis- 
order with quantum phase fluctuations plays a central 
role for the superconductor normal-metal phase transi- 
tions in this class of quasi-lD superconductors. Recent 
experimental data found in Refs.9,10 are shown to be 
consistent with the present theory. Quantum criticality 
is controlled by two quantities, namely disorder strength 
and a dynamical parameter of phase fluctuations. The 
present model's quantum criticality signals the existence 
of a quantum critical phase between SC- and normal 
phase. Its nature deserves further investigation. 

In our study we neglect the effects of non-linear ex- 
citations, which are a subject of current interest in low 
dimensional systems. Note that this topic was explicitly 
studied by us for quasi-2D Josephson coupled supercon- 
ductors in Ref. 38. As we have shown in the previous 
section, the classical motion of the phase is described by 
a system of coupled sine-Gordon type non-linear equa- 
tions (24) , which contains non-linear dynamic excitations 
as well as static topological defects. The self-consistent 
phonon approximation allows us to calculate the phase- 
phase correlator between two arbitrary points r = {z,}} 
and r' = {z',j + g} 

((cos(</.j(^) ^ = e-^«(^-^''^'^), 

which can be shown to decrease at g — >■ and \z — z'\ — )■ 
oo as a power law ^ (^n/l^ ~ z'\)^^'^\ setting up a 
quasi-long-range order and implying the existence of a 
Berezinskii-Kosterlitz-Thouless (BKT) topological phase 
transition''^ (perhaps at T = 0) in a single SC wire. The 
critical index f3(T) contains both phonon and vortex con- 
tributions. Although the phonon contribution to /3{T) 
can be calculated within the SCHA, it is not clear how 
the vortex contribution changes the former one. In our 
knowledge, the mechanism of excitations of the vortices 
with opposite fugacities and their binding in quasi-lD su- 
perconductors has not been adequately studied, and the 
topic needs further investigations. 
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VII. APPENDIX 



In order to calculate Sa{g,0) for g = BxN^ + GyNy 
(where e^jSy are unit vectors and Nx,Ny are the num- 
ber of unit cells in directions of x, y, correspondingly) we 
rewrite Eq.(34) in the following form: 

c / n^ '^1^ r '^1^ r '^ly 1 ^ cos(q_Lg) 

7-1 ^71^ J-TT 27r J_^ Ztt . /„2 j_ X2 ('n^,.,2 



(56) 

where uj'^{qx,qy) = 2(2 — cosg^ — cosQy) and qj^ = 
{Qx,Qy}- Introducing a new variable z = qz/Squ{0) and 
using the transformation 



1 _2_ r 



■ 2 / 2 I 2 \ 



(57) 



one can integrate out z,qx and qy in Eq.(56). Finally, 
5c (g, 0) is expressed as an integral over u = 2t^ as 



5a(g,0) 



2 Jo u' 



xerf(5-J(0)VV2), (58) 



where In{u) is the Bessel function of an imaginary ar- 
gument, and erf (2:) = dte~* is the error function. 
By using the following asymptotic expressions for In{u), 



I.{z) 



Eoo 
k=0 



fc!r(y+fe-|-l) \2) 

(^-H/2)(./-l/2) j 



and for erf (2;) 



erf(^) 



< 2; < \/z7TT 
z»v 



z < 1 



(2fc-l)(fe-l)! ' 
fe(2fe-l)!! ^ 



we get for, e.g. Sa.{N, 0) at = N,Ny = the following 
expression 



Sc{N,0) = aln 



Ci 



5,«(0) 



a C2 
4^ iV' 



(59) 



i.e. e-^^"W") ^ (<5,„(0))"exp(||^), where Ci and 
C2 are constants of order of unity. If we take only the 
first terms in the expansions of In{z) and erf (2), we get 
Ci = 1 / \/2 and C2 = 1 . Higher order contributions cor- 
rect only these constants. Thus it is seen from Eq.(59) 
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that the phase-phase correlator in the transverse direc- 
tion saturates at T = to its asymptotic value of S"^{0) 
for distances of the unit cell size a± . 

The correlator Sa{g,T) at T 7^ is also calculated in 
the (g, 0) was obtained above. Using in 

Eq.(40) the representation coth(7rz) = f J2'^=_^ n'^+z'^ ' 
the correlator Sa (g, T) assumes the following form 



Sa{g,T) = 



(1) 



E 



dqz 



27r 



dqy 



1 - cos(q_Lg) 



2tt 



(60) 



We use the transformation 1/a = e~°-^du in Eq.(60) 
and carry out the integrations over qz . q.,- and qy . For the 
particular case of = N and Ny = 0, the expression 
for Sa{g, T) is reduced into the following form 

S„{N,T) = ^ r ^e-2«/o(u)[/o(«) - iNiu)] X 
W'^ Jo V" 



cerf [6-,\T)^)<l>lu, 



(1) 



(61) 



where $(u, r) is given by 



$(m, r) = r exp{— ur^n^} 

n=—oo 

/ T(l + 2e-"^'), u>T- 

U <T~ 



(62) 



here r represents the normalized temperature (see 
Eq.(54)) and the sum in Eq.(62) is also known as the 
so-called EUipticTheta-function ^3(0, exp(— mt^)),"*^. Us- 
ing the asymptotic expressions for the Bessel and the 
error function as well as for $(u,t) in Eq.(61), we get 
the following explicit expressions for exp{—Sa{N,T)} 



«?„exp{-C^ + -ife^}. < 25^ < 1 < < JV 



where C3, C4, and C5 are again constants of order unity. 



^ D. Jerome, A. Mazaud, M. Ribault, and K. Bechgaard, J. 

Physique Lett. 41, L95 (1980). 
^ C. Bourbonnais and D. Jerome, arXiv:0904.0617vl [cond- 

mat.str.el] 3 Apr. 2009. 
3 M. -Y. Choi, P. M. Chaikin, S. Z. Huang, P. Haen,E. M. 

Engler, and R. L. Greene, Phys. Rev. B 25, 6208 (1982). 
* S. Bouffard,M. Ribault, R. Brusetti, D. Jerome, and K. 

Bechgaard, J. Phys. C 15, 2951 (1982). 
^ B. J. Powell, J. Phys.: Condens. Matter 20, 345234 (2008). 
^ A. A. Abrikosov, J. Low Temp. Phys. 53, 359 (1983). 

J. Shinagawa, Y. Kurosaki, F. Zhang et al., Phys. Rev. 

Lett. 98, 147002 (2007). 
® T. Ishiguro, K. Yamaji, and G. Saito, Organic Supercon- 
ductors (Springer- Verlag, Heidelberg, 1998). 
® N. Joo, P. A.-Senzier, C. Pasquier, P. Monod, D. Jerome, 

and K. Bechgaard, J. Euro. Phys. B 40, 43 (2004). 
° N. Joo, P. Auban-Senzier, C. R. Pasquier, D. Jerome, and 



K. Bechgaard, Europhys. Lett. 72, 645 (2005). 
" B. J. Powell and R. H. McKenzie, Phys. Rev. B 69, 024519 
(2004). 

1^ P. W. Anderson, J. Phys. Chcm. Solids 11, 26 (1959). 
A. A. Abrikosov and L. P. Gor'kov, Zh.Eksp.Teor.Fiz. 39, 
1781 (1960) [Sov.Phys.JETP 12, 1243 (1961)]. 
V. M. Galitski and A. I. Larkin, Phys. Rev. B 66, 064526 
(2002). 

1^ C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000). 

1^ P. Lee, Phys.Rev. Lett. 71, 1887, (1993). 
^'^ P. J. Hirschfeld and N. Goldenfeld, Phys. Rev. B 48, 4219 
(1993). 

^® Y. Suzumura and H. J. Schulz, Phys. Rev. B 39, 11398 
(1989). 

^® E. P. Nakhmedov and E. V. Tahirov, J. Phys.: Cond. mat- 
ter 6, 2245 (1994); ibid 6 2259 (1994). 



12 



^° E. P. Nakhmedov and M. Hiiner, Physica C 299, 301 

(1998) 

21 M. Ma and P. A. Lcc, Phys. Rev. B 32, 5658 (1985). 

22 R. Oppermann, Physica A 167, 301 (1990). 

23 A. Ghosal, M. Randeria, and N. Trivedi, Phys. Rev. B 65, 
014501 (2001). 

2" B. Spivak, A. Zyuzin, and M. Hruska, Phys. Rev. B 64, 
132502 (2001). 

25 V. M. Gahtski and A. I. Larkin, Phys. Rev. Lett. 87, 
087001 (2001). 

2"^ S. Sachdev, P. Werner, and M. Troyer, Phys. Rev. Lett. 

92, 237003 (2004). 
2'' B. Spivak, P. Oreto, and S. A. Kivelson, Phys. Rev. B 77, 

214523 (2008). 

2® A. D. Maestro, B. Rosenow, M. Miiller, and S. Sachdev, 

Phys. Rev. Lett. 101, 035701 (2008). 
2® F. Mondaini, T. Paiva, R. dos Santos, and R. T. Scallettar, 

Phys. Rev. B 78, 174519 (2008) 
3° R. Oppermann, Z. Phys. B - Cond. Mat., 63, 33 (1986) 

31 A. M. Finkel'stein, Pis'ma Zh. Eksp. Teor. Fiz. 45, 37 
(1987) [JETP Lett. 45, 46 (1987)]. 

32 P. W. Anderson, K. A. Muttalib, and T. V. Ramakrishnan, 
Phys. Rev. B 28, 117 (1983). 

33 M. A. Skvortsov and M. V. Feigel'man, Phys. Rev. Lett. 
95, 057002 (2005). 

3" K. B. Efetov and A. I. Larkin, Zh.Eksp.Teor.Fiz. 66, 2290 

(1974) [Sov.Phys.JETP 39, 1129 (1975)]. 
35 Yu. A. Firsov and G. Yu. Yashin, Zh.Eksp.Teor.Fiz. 72, 

1450 (1977) [Sov.Phys.JETP 45, 761 (1977)]. 
3« H. J. Schulz and C. Bourbonnais, Phys. Rev. B 27, 5856 

(1983). 



3'' V. J. Emery and S. A. Kivelson, Nature 374, 434 (1995); 

Phys. Rev. Lett. 74, 3253 (1995). 
3« E. P. Nakhmedov and Yu. A. Firsov, Physica C 295, 150 

(1998). 

39 T. M. Rice, Phys. Rev.140, A1889 (1965). 
^° The ground state of a non-interacting electron gas in an iso- 
lated ID wire or in a quasi-lD system with a nested Fermi 

surface is a Peierls dielectric. Suggestion of n- subband in 
the cross-section of a pure single wire surves to prevent the 
Peierls transition in the system. Disorder in the quasi-lD 
system destroys also the perfect nesting, suppressing at the 
same time the Peierls transition. If isolated wires contain 
point-like impurities as well, the ground state becomes a 
dielectric phase due to Anderson localization, and the tran- 
sition becomes superconducting-insulator phase transition. 
In our case, anion disorder in the organic superconductors 
may not modify the Fermi surface. 

"1 According to Eq.(33) |Aq|2 oc (1/iV), this approxima- 
tion is valid for N ^ 1. Strictly speaking, this har- 
monic approximation is valid for small temperatures when 
Sq l^q|'^{-^q) < 1) which means that few phonons are ex- 
cited in the system. 

''2 V. L. Pokrovskii and G. V. Uimin, Zh. Eksp. Teor. Fiz. 65, 
1691 (1973) [Sov. Phys. JETP 38, 847 (1973)]. 

'*3 A. I. Larkin, Pis'ma v Zh. Eksp.Teor.Fiz. 2, 
(1965) [Sov. Phys. JETP Lett.2, 130 (1965)]. 
Q. Yuan, H.-Y. Chen, H. Won, S. Lee, K 



205 



and C. 



Maki, P. 

S. Ting, Phys. Rev. B 68, 174510 



Thalmeier 
(2003). 

See, e.g. P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987) 
Mathematica handbook. 



